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Abstract. We define a discrete gauge-invariant Yang-Mills-Higgs action on spacetime simpli- 
cial meshes. The formulation is a generalization of classical lattice gauge theory, and we prove 
consistency of the action in the sense of approximation theory. In addition, we perform numer- 
ical tests of convergence towards exact continuum results for several choices of gauge fields in 
pure gauge theory. 
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1. Introduction 

This article is the natural prolongation of an earlier work [2^ . In [2j we defined a gauge invariant 
discrete action to approximate the Yang-Mills-Higgs (YMH) action |17| [T3l [M] [T2] on a simplicial 
mesh approximation of a spatial domain, and we proved consistency of the discrete action for 
smooth fields. In this article we expand the discussion to 4d spacetime simplexes in time and 
three spatial dimensions, and we also expand the consistency proof to fields in the natural energy 
norm. 

More precisely, we define the spatial part of the discrete action as in [2, , using Whitney forms 
|15j . Wilson lines and loops |T6|l6|. and concepts from FEM [5l|9j[10]. We then expand the spatial 
mesh to spacetime by uniformly discretizing time. The spacetime mesh is defined by repeating 
the spatial mesh at every time step. Furthermore, we extend the Whitney forms to spacetime, 
yielding a natural expansion of the discrete spatial action to a spacetime action. 

The purpose of this work is to prove consistency of the resulting approximation scheme for the 
classical YMH equations, and also to develop an analogue of lattice gauge theory, i.e. performing 
quantum field calculations using the introduced action. The latter will be described in a companion 
article |8]. 

Key words and phrases. Lattice gauge theory, QCD, simplicial complex, Yang-Mills theory. Finite element 
method, MSC2010: 35Q40, 65M50, 74S05, 81T13, 81T25. 
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2 TORE GUNNAR HALVORSEN AND TORQUIL MACDONALD S0RENSEN 

We use the same notation as in [2 , and the article is organized as follows. In section [2] we 
introduce the continuous YMH action, and also develop the proposed siniplicial gauge theory. In 
section [3] we review some properties of the differential of the exponential map for matrix groups 
and the Baker-Campbell-Hausdorff formula. We prove consistency of the action in the energy 
norm in section |4] Finally, in section [5] we describe some numerical convergence tests, performed 
on several gauge fields for which the continuum action was exactly calculable. 

2. Construction 

2.1. The continuous Yang-Mills action. Let M = M x 5 be a riemannian spacetime manifold, 
where M represents time and 5 is a bounded domain in three dimensional euclidean space. We let 
M be equipped with a lorentzian or euclidean signature and coordinates x = (t,x). Furthermore, 
let Q he 'A compact Lie group with associated Lie algebra g, and assume that Q can be represented 
by a subgroup of the complex unitary n x n matrices, for some n. The hermitian conjugate of a 
matrix g is denoted g^ and the real valued scalar product is 

(2.1) g'-g:=mrig'g"). 

The space of smooth fc- forms on M is denoted ri'^(M), and the space r2'^(M)(g)g can be identified 
with the space of smooth g- valued fc-forms on M. The bracket of Lie algebra valued forms is defined 
as 

(2.2) [u^g,u'®g']:={uAu')® [g, g'] , 

where u, u' are real valued differential forms and g, g' E g. 

A connection one-form on M is an element A — {Aq,A) G il^{M) (g) g, where Aq represents 
the time component and A the spatial component. The temporal curvature J^*(A) and spatial 
curvature T^ (A) of such a one- form are given by 

(2.3) J'\A)^dAo + dtA + [Ao,A], J^'iA) ^ dA + hA,A], 

where dt and d denote exterior derivative in the temporal and spatial direction respectively. We 
will also use such forms with lower regularity, see e.g. [3i. The independent variable in Yang-Mills 
theory is such a connection one-form (not necessarily smooth) , and the action describing it is given 

by S[A] = S[A]t + S[A],, where 

(2.4) S[A]t := f |J-*(A)^ 5[A], := / \^'{A)f. 



A gauge transformation of a connection one-form A is associated with a choice of G G ^ for 
each a; € M, such that 

(2.5) Ao{x)^G(x){Ao{x) + dt)G-'^{x), A{x) ^ G{x) {A{x) + d)G~\x). 

The action S[A] is invariant under such gauge transformations. 

2.2. The interpolated FEM action. In the FEM formulation, we assume T to be a siniplicial 
complex spanning the spatial domain S. The simplexes are referred to as vertexes, edges, faces 
and tetrahedra according to dimension, and are labeled «, e, / and T respectively. The symbol T 
will also be used for simplexes of any dimension. We also suppose an orientation has been chosen 
for each simplex in T. In addition, we assume that time is discretized with a time step At, and 
that the siniplicial complex T is repeated at every time step, resulting in a spacetime siniplicial 
complex T. For a more detailed definition of simplicial complexes, consult [4, section 5]. 

Thus, as the basic building block in classical FEM theory is a tetrahedron T, the basic building 
block in this extended FEM version is T x /,-, where Ir = [t,t + At] and r denotes the temporal 
nodes. 

Furthermore, let W''{T) {W''{T)) be the space of Whitney k-forms [15j on T (T), with canonical 
basis (Ay), T ranging over the set of /c-dimensional simplexes in T. The 0-forms A^ are the 
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barycentric coordinate maps taking the value 1 at vertex i and at all others. For an edge 
e = {hj}, the associated Whitney 1-forni is defined by 

(2.6) Ae — XidXj — XjdXi, 

and for a face / = {i,j, k} the associated Whitney 2-form is defined by 

(2.7) -^z = 2 {XidXj A dXk + X^dXk A dXi + XkdXi A dXj) . 

In order to formulate the Yang-Mills theory in a spacetime FEM setting we need to extend these 
k-forms to k-forms on T. In addition, we need to define the temporal edge and temporal face basis 
functions, which are constructed as in |3]. 

The spatial Whitney k-forms are extended to be piecewise affine in time and are denoted (Ax(t))i 
i.e. 

(2.8) At^ AT(r) = At«)Pi*, 

where P^ denote polynomials in the time variable of degree at most one, and T{t) := (r, T) denotes 
the spatial simplex T at temporal node r. More precisely, Ay(^) is the piecewise affine function 
in time, taking the value At at r and at all other temporal nodes. This is consistent with the 
requirement that the tangential part of the curvature of the gauge potential should be continuous 
across faces. 

The temporal edge basis functions are constructed as follows. To every vertex i in the spatial 
mesh, there are temporal edges et(T) — {VjV-i-At}, where ir '■— i{T). The temporal basis edge 
function attached to et(T) is then the piecewise constant function in time defined by 

(2.9) A.^(.)(t) = ^'°"^'''' '^'^ 

10, otherwise. 

Here, tt is the canonical projection onto S, i.e. 

(2.10) tt:M = M.-x.S^ S, 

and dt is the standard basis one-form in the temporal direction. 

Finally, the temporal face elements are constructed as follows. To every spatial edge e there 
are corresponding temporal faces /((r) — e x Ir- Consider the spatial Whitney edge element Ag. 
Then apply the pull-back of tt to construct a one-form on spacetime, and then wedge it with dt. 
More precisely, denote by 

(2.11) TT* : n{s) -^ n{M), 

the pull-back induced by n. Then the temporal basis face function is 

(2.12) U*iX^)Ai-^dt^X^o.Ai-^dt, tel^ 

^ ' ^*^ -"^ ^ [O, otherwise. 

This construction ensures that the temporal face basis is orthogonal to the spatial face basis. 
The space spanned by (Ai^r)), {K{t)) ■= {iKtiT),K{T))) and (Af(^)) := {{Af^(^r),-^f{T))) are 
denoted W''(T), W^(T) and W^(T) respectively. If no confusion can arise, the time dependence 
index r that identifies a temporal node is usually omitted to compactify notation. 

Thus, let A = {Ag, A) G W^ (^g, A — J^e ^o,et Ae^ + J^e ^eAg, where the summations are over 
oriented edges, and we remark that 

(2.13) Ao^e,^ f Ao, A,= [ A. 

J et J e 

The temporal and spatial curvatures of A are given by 

J"* (A) =J2^edtAe+J2^0,etdAe, + ^ [Ao,e, , ^e] Ag, A Ae, 
e et et,e 

(2.14) , 

T'{A) = ^ AedAe -f -^K,^e']Ae A Ae'. 
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Since A^ AA^' ^ W^(T) we choose to interpolate J"* (A) and T^iA) onto W^(T), instead of working 
with higher order Whitney elements. 

Let (/*,/*) and (J*, J'*) denote interpolation operators onto the temporal and spatial Whitney 
one- and two- forms, respectively. They are projection operators defined by 

(2.15) 

ft ''f^ f •'^ 

and are well defined in particular as maps ^^(M) -^ Wi(T) and n'^{M) -^ W2(T), respectively. If 
we define /:=/*+ 1'^, J := J* + J*, and d := {df, d), we have do/ = Jodby Stokes' theorem. 
In particular, do P ~ J^ o d, consistent with the classical Whitney elements. 
The degrees of freedom of the interpolated curvatures are 

(2.16) 

4(J-*(A)) := / F\k) = ^ ^e + 51 ^0,«, + Y. Ce,e[^O.e„ Ae], Ce,e = / Ae, A A^, 

•^' eeft eteft e,eteft ^* 

J]iT^A)):= f ^%A)=J2Ae+l J2 Cee'[Ae,A,,], C,,. ^ f A, A A,, . 



eef e.e'ef 

The interpolated FEM action is then S-^[A] = S-^[A]t + S-'[A]s, where 



t = / I J^-F'(A)|^ = 

Ju 

- 5i E M/,/,tr (4(^*(A))j;,(^*(A))«) , Mf^f, := f A,, • A,,, 

(2.17) 

S''[Al = / |J^J-^(A)p = 

Ju 

= 5R^ A%,tr {JI{r^{A))jp{:F%A))") , Mjf, := [ Af ■ Af,, 

f,f' -^M 

where (•) denotes the scalar product of alternating forms w.r.t. the Lorentzian signature. 

By inspiration from lattice gauge theory |16l E], we will now in several steps construct an 



approximation to this action that is gauge invariant, i.e. invariant under the transformation (2.5 1. 

2.3. An intermediate action. 

The spatial part. Let {i,j,k,l} be the vertices of T, and pick A € W^(T) (E) g. Attached to an 
edge e = {i,j} oriented from i to j one has an element Ae = Aij e g, where we recall that 
Ae — J A, and parallel transport from i to j is given by C/^- = exp(Ajj). We suppose Uij to be 
close enough to 1 so that its logarithm is unambiguous. We use the sign convention Aij — —Aji, 
which corresponds to Uji = U~j^. 

The discrete spatial curvature associated with a face / = {i,j, k} is defined by 

(2.18) i^% := U.,,U,,Uk^. 

and is in analogy with square Wilson loops in classical lattice gauge theory [16]. By definition, 
this formula locates the curvature at vertex i. The curvature at vertex j is related to this by 

(2.19) ^, = C/,,i^,[A„ 

which gives a formula for parallel transport of curvature from i to j. Concerning orientation of a 
face, we notice 

(2.20) ^% = (F4^-)"^ 



SIMPLICIAL GAUGE THEORY ON SPACETIME 5 

When / is oriented i ^^ j ^^ k, and the curvature is located at i, we write Ff- — -F/.j.. Thus, 
we have defined the spatial curvature of a pointed oriented face /. The distinguished point of / is 
denoted /, which in this example is f — i. 

A discrete gauge transformation is associated with a choice of Gi Cz G for each vertex i. One 
then transforms A such that 

(2.21) U,j^G^U^,G-\ 
implying that the spatial curvature transforms as 

(2.22) F^ ^ G,F^G-\ 

similarly to the gauge transformation of the field strength tensor in continuous Yang-Mills gauge 
theory. 

The temporal part. Concerning the temporal part of the curvature, the construction is similar. Let 
{(ir, jTiir+At, ir+At} be the vertices of the temporal face ft{T), and pick A ~ {A^, A) e W"'^ (E) g. 
Attached to a spatial edge e(T) = {i^^j^.} one has as before an element Ag^^-j —: Ai^j^ G Q, and 
parallel transport from v to jr is again given by Uijlr) = exp(y4.i^j^). 

Attached to a temporal edge et(r) ~ {irjir+At} one has an element Ao,ej(r) = ^o,i^i^+At ^ 0; 
where we recall ^o,et = / Aq, and parallel transport from i^ to ir+At is given by C/o,i^i^^^j = 
exp(^o,iTiT+At)- ^^ ^^^ ^^^ sign convention ^o.JTix+Af ~ ^^Q-iT+^ti-r which corresponds to 
C/o.i^iT+At "^ ^(Ti i ■ Again, we suppose Uij and C/o.i^i^+^t to be close enough to the iden- 
tity, so that their logarithms are unambiguous. 

The discrete temporal curvature associated to /((t) is 

(2-23) ^l,,.,,.+A.v+A. = U,,{t)U,,,^,^^^,U,,{t + Ai)C/o,v+A.v, 

and again we see the similarity with classical lattice gauge theory. This formula locates the 
temporal curvature at vertex v. The curvature at vertex ir+At is 

which gives a formula for parallel transport of curvature from i^ to ir+At- Concerning the orien- 
tation of a temporal face, we notice 

^ -' iT.JT.ix + At.ix + At ~ ^ ir , ir + At Jt + At Jt'' 

When /((t) is oriented v — > jr — > jr+At -^ ir+At, and the curvature is located at v, we write 
Fl ,. = Fi^j^j^^^t,i^^^t- The distinguished point of this pointed oriented face /((r) is denoted 

Mr). 

Under a discrete gauge transformation, one transforms Aq such that 

(2-26) ^o,irir+Af ^ Gi^UoA^i^^^tG^^^^^, 

implying that the temporal curvature transforms as 

(2-27) f;^(,)K^G.^F;^(,)G71. 

The action. We define the intermediate action as S^[A] — S^[A]t + S^[A]s, where 

(2.28) S'[A]t = 5ft ^ M/,/^,tr [(^;.^(A) - 1)(^;,(A) - 1)" 

ftJi 

and 

(2.29) S'[Al = 5ft^ M^-;,tr [{Ff{A) l)(Ff,{A) 1)"] 

fj' 
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2.4. The simplicial gauge theory action. In this section we use the same notation as in the 
construction of the intermediate action, equations (2.281 and (2.291. 

The goal is of course to formulate a gauge invariant action. The intermediate action is not gauge 
invariant due to interactions between different faces with non-coincident distinguished points. 
However, this can be resolved by parallel transport. We treat the spatial and temporal part 
separately. 

The spatial part. Let / and /' be two spatial faces of a tetrahedron T. The spatial curvature 
at / and /' can be connected by the parallel transport operator Ufr,. However, the curvature 
associated to the face / at time t will interact with the curvature associated to the face /' not only 
at time r, but also at times r ± At. This is resolved by the parallel operator in time U^ i,^-.i,,-.. 
Thus the troubled term 

(2.30) (F;(,) - l)(^f (.,) - 1)^ 
is replaced by 

(2.31) C/^,(^)y(^)(F;(^) - l)t/y(^)y,(^)f/oy,(^)y,(^,)(-F;,(^,) - 1) C/q y,(^,)y,(^), 
and from the transformation properties of F and U , its trace is now gauge invariant. 

The temporal part. Let /((r) and f[{T') be two temporal faces. By properties of the canonical 
basis (AyJ we know that interactions between the temporal curvature occur only at equal time 
intervals. Thus the term 



^FUr)-^){Fl^(r-)-^f 



(2.32) 

is replaced by 

(2.33) [/y,(^)y^(^)(F;^(^) - l)C/y^(^)y,(^)[/oy,(^)y,(^,)(F;.,(^,) - l)"U^J,^^^,y,(^^y 

and we note that r = r' or r = t' + Ai. By the transformation properties of F and U , the trace 
of this term is also gauge invariant. 

The action. We define the simplicial gauge theory (SGT) action as S^[h] — S^[k\t + S^[A\s where 



s' 



^Mr) 



^Mrjiiir)^ 



(2.34) 



and 



(2.35) 



X U, 



OJi(r)mr') 



^hir') 



H 



u, 



aJl{r')n{r)J' 



/M,/'(r') 

Xt/n 



^/V) " 



^/(r)/'(r)X 



^kr') 



1 



H 



U, 



aj'{T'}f'(T) 



^0,/'(r)/'(r') 

Thus we can conclude. 

Theorem 1. The simplicial gauge theory action S^ is discretely gauge invariant. 

2.5. Gauge invariant scalar field action. To complete the Yang-Mills-Higgs-action, we add a 
discrete gauge invariant action for the scalar Higgs field (however no consistence proof for it will 
be provided) . The gauge group Q can be represented by a subgroup of the complex unitary n x n 
matrices, and in that representation the basic scalar fields will form an n-tuplet, i.e. 

/ M^) \ 

(2.36) X K> (j){x) = : e C". 

\ K{x) J 
The action describing this multi-component field is given by ^[i?!). A] — S[(f>,AQ\t + S[(l),A]s, where 



(2.37) 



S[cl),Ao]t^ 



I \DaA\ S[c^,A]s= [ \DA<f>\', 

JM Jm 
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and DjigCJ) — {dt + ^o)0 ^iid Da(J) — {d + A)(p are the gauge-covariant derivatives on scalar fields. 
We observe that the action is invariant under the following set of local gauge-transformations 

(j){x) n- G{x)(j)(x), 

(2.38) Ao{x) ^ Gix) (Aoix) + dt) G~^{x), 

A{x) ^ G{x) {A{x) + d) G-\x), 

where G{x) G G- This is so, since the covariant derivatives transform as Dao4>{x) i->- G{x)DAo(f>{x) 
and DA(l>ix) ^ G{x)DA(i){x), and G"G = 1. 

Our aim is to construct a FEM inspired action for this scalar field that is gauge invariant, and 
the key ingredient is again the parallel transport operator U introduced in previous sections. 

The scalar field is formally a zero-form, and in FEM it has degrees of freedom at the nodes of 
the mesh. Thus, let G W" (g) C". Then we can write 

(2.39) = ^0,^A,^, </,,^=0(v). 

The covariant derivatives of cj) are one-forms, they have degrees-of-freedom on the edges of the 
mesh, and are approximated as in lattice gauge theory (LGT) |16| f6|. More precisely, let e* be 
a temporal edge and e a spatial edge of the mesh. Their origins are denoted e-t and e, and their 
targets et and e, respectively. Then the components of the covariant derivative of (f> along these 
edges are approximated as 



(2.40) 



We observe that 

(2.41) {5A,4>)e,^Ge,{5Ao^)e„ {5 a4>) e ^ G ^{5 A^) e 

whenever 

(2.42) ^,^ ^ G,^cb,^, f/o..+.... ^ G,^+,,{/o.^+,,.^G-\ {/,, ^ G^U.^GJ^ 

Thus the components of the approximated covariant derivatives transform as in the continuous 
case. However, this is not enough to ensure local gauge invariance of the action, since the inner 
product of edge basis functions involves interactions between different edges. 

The temporal part. By using the approximation of the temporal covariant derivative in equation 



(2.401, a FEM inspired approximation of S'[(^, Ao]t reads 



(2.43) S^[<l>,Ao]t = 3ff V Me,e;(5Ao0)f!('5Ao0)ei, Me.eJ = / Ae, • A, 



We note that this approximation is not gauge invariant since the mass matrix M^^e' is not diagonal. 
This problem can be overcome by either mass-lumping of Mf.^e' lUi or by using the parallel trans- 
port operator U to localize the terms in the action. The mass-lumping procedure severely restricts 
the structure of the mesh for time-dependent problems, thus we choose the second alternative. 
We make the replacement 

(2.44) {SAoCl^)^MAo<t>)e', ^ {SAo<l>)ZUe,e'MA,cl,),'^, 



in equation (2.431, and approximate the temporal part of the action as 



(2.45) 5^[0,Ao]t=3fi J2 M,^{r)e[(r){6AA)"ir)U,,{r)e',(r){SAA)e',{r)- 

By the transformation properties of (Sao^') ^nd [/, this is gauge-invariant. 



TORE GUNNAR HALVORSEN AND TORQUIL MACDONALD S0RENSEN 



The spatial part. By using the approximation of the spatial covariant derivative in equation ( 2.40 1 , 
a FEM inspired approximation of S[(l), A\s reads 



(2.46) 



5^[0,^], =5R^M,e'(<5A'^)f(<5y 



M,. 



K ■ Ap/. 



Again, this approximation is not gauge invariant since the mass matrix M^e' is not diagonal, and 
again we choose to resolve this problem by parallel transport. 
We make the replacement 

(2.47) 



{^A(t))e(r){^A(t))e'(T') ~> ('5A'/>)e(r) ^0,e(r)e(r') ^e(r')e'(T') (S A(t)) e' (r' ) 



in equation (2.461, and approximate the spatial part of the action as 

(2.48) S^[(l),A]s^di ^ Me{r)e'{r'){SA4')e{T)Uo,S{r)e(T')Ue(r')e'{T'){SA<l))e'{r')- 

e(T),e'(r') 

By the transformation properties of {6A(t>) and U, this is gauge invariant. We therefore have 
Theorem 2. The action S^[(f),A] = 5^[0, Aojt + S^[(j>,A]s is discretely gauge invariant. 

3. The differential of the exponential map for matrix groups 

Let Q he a, compact matrix Lie group with associated Lie algebra g = ^Td^- '^^^ structures of 
G (the connected component containing the identity) and g are related through the exponential 
map 

(3.1) exp:g^g, ^ i-^ exp(A), 

which for matrix Lie groups is given by the usual power series expansion 



(3.2) 



exp(A) = Y. 



A" 



In order to apply Hamilton's variational principle to the gauge invariant SGT action introduced 
above, we need to calculate the differential of the exponential map in an arbitrary direction. In 
other words, we need to find an expression for 

(3.3) i^(e^) • B '^ 



^A+tB 



dr 



r=0 



If A and B commute, i.e. [A, B] := AB — BA = 0, the differential is straightforward to calculate 
and equals 

(3.4) D(e^) • B = e^B. 

This is not the case if [A, B] ^ 0. 

However, one can prove the following useful formula [T 

Proposition 3. Let X and Y be n x n (n £ N) complex matrices. Then the following relation 
holds 

- adx 



(3.5) 



o-^+T-y 



dT 



X 



T = 



1 



ad 



Y 



X 



where adx ^ ~^ OLi^d adxY = [X,Y]. More generally, if t i— > X{t) is a smooth matrix evaluated 
function, then 



(3.6) 



dr 



^X(r) ^ ^Xir) 



I dr 



ad 



'X{t 



By the cyclic invariance of the trace operator, i.e. if A, B, C are nxn matrices then tT{ABC) = 
tr{BCA) = ti-{CAB), one can also prove 

d 



(3.7) 



dr 



tr(e 



X+tY-' 



T = 



tr(e^y). 
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3.1. The Baker-Campbell- Hausdorff formula. For complex numbers x and y we know that 

(3.8) e^e^ = e^+^. 

This is not the case when x and y are replaced by matrices X and Y , but the following relation 
holds Umi]: 

Proposition 4. Let X and Y be n x n (n £ N) complex matrices. Then 

oo 

e^e^^e^, Z==^c„, ci=X + Y, 

(3.9) ["/2] 

(n+l)c„+i = -[X-y,c„]+^-^ ^ [cfe,,[...[cfc,^,X + r]...]], 

fci+...+fc2p=n 

where [n/2] denotes the largest integer snaller than n/2, and Bj is the j-th Bernoulli number. 
The first four terms in the expansion of Z is given by 

ci=X + Y, C2 = l{XY-YX), 

(3.10) C3 = —{X'^Y + YX^ - 2XYX + Y^X + XY^ - 2YXY) , 
C4 = TT^i^^Y^ - 2XYXY - Y^X^ + 2YXYX) . 

4. Consistency 

Since we have formulated the theory in a spacetime FEM setting, we choose to define consistency 
for the entire action, not only for the spatial part as is usual. This definition of course encompasses 
the usual one. Note that we only prove consistency for the gauge field action, without the scalar 
field. Inclusion of the scalar field is a simple extension of this proof. 

We suppose that we have a regular sequence of simplicial meshes T^ of the spatial domain S. 
The diameter of a simplex T is denoted Ht, and the biggest Ht when T e 7^ is denoted ft.„. In 
addition, we suppose that time is discretized by a time step Ai„, and that Tn is repeated at every 
time step, resulting in a simplicial mesh T„ of the spacetime domain M. We suppose that 

(4.1) max{(/i„),(Ai„)}"^0. 

The interpolation operators onto the Whitney elements introduced earlier are attached with a 
subscript n to associate them with the mesh T„. Finally, let X„ = W"'^(T„) (g) g. 

Definition 5. We say that two actions Sn and S'^ defined on X„ are consistent with each other, 
with respect to a norm || • ||, if for any A G Xn we have 

\DS,AA]A' ^ DS'JA]A'\ 
y ' °^t^ \\AHi ^T^ 

If there is a constant C not depending on n such that quantities a„ and &„ satisfy a„ < C6„ 
for all n, we write a„ = 0{bn). To compactify notation the subscript n will be suppressed. 
We have introduced three different actions S''^[A], 5'^ [A] and S''^[A], and the plan is to show 

(1) 5'"' [A] consistent with S[A], 

(2) ^^[A] consistent with S-^[A], 

(3) S^[A] consistent with S^[A], 

which implies the consistency between iS'^[A] and ^[A]. 
We will prove consistency in the energy norm, i.e. 

(4.3) Ao,AeL°-{R;H\S)), 9*^0,9*^ e ^""(M;^^ (5)), 
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where dt is a shorthand for d/dt. To compactify notation we define 

■ \\lp(R-Li{S))^ VO <p,q < oo, 



(4.4) 



II ■ \\LP{Li) 
■ llL-(ifi) 

l|A|| 



l^ol 



f;ffl(S))) 

l^ll- 



The H^ spacetime EucHdean seniinorni is denoted | • |//i(m)- 

4.1. Consistency between 5"' [A] and ^[A]. The two actions are given in equations (2.17) and 

(2.4) respectively, and 

\DS[A] ■ A' - DS'^[A] ■ A'l < Kt + K„ 

(4.5) Kt := \DS[A]t ■ A' - DS-'[A]t ■ A'|, 

Ks := \DS[A]s ■ A' - DS'^[A]s ■ A'\. 

We treat Kt and Kg separately. To compactify notation we define 

cIaA' := dA'o + dtA' + [A'o.A] + [Ao,A% 



(4.6) 



dAA' 

dA 

dAo 



^dA' + [A,A% 
= dtA + WA, 
= dtAo + WAo, 



where V is the spatial gradient on k- forms. 



(4.7) 



We first estimate the term Kt. From equations ( |2.17| and ( |2.4| we get 
Kt^ 



((J"*(A),dAA') - (J*J-*(A), JMaA') 

< 2||^*(A)|U2(i2)||dAA' - J*dAA'|U2(i2) + 2||^*(A) - J*J-*(A)|U2(i2)|| J*dAA'|U2(i2). 

The interpolation operator J* possesses two important properties. First of all it is a projection 
operator, and second it is stable L^ — >■ L^ by scaling. This implies that 

\\<d^A'-j'<d^A'\\L2(^L^)<Ch{\[AlA]\„i(M) + \[Ao:A']\„i^M)) 

< Ch {\\A',dAU2(L2^ + \\AdA'„U2^L-) + WA'dAoWmL^ + WAodA'h^^L^) 

(4.8) + P'||L~(L6)l|aAo||L»(L3) + ||Ao||L»(L6)||a.l'|U^(i3)) 

< C/ll/2(P^|Uoo(i.)||9A|Uoo(i2) + \\A\\Lo.^Le)\\dA'o\\L^^L2) + 

+ ||A'||ioo(i6)||aylo||L~(L2) + ||Ao||L-(L«)||aA'||ioc(i2)) 

< C/ii/2 (||A'|U=.(^i) + ||9,A'|U.o(i2)) 

where C is a generic constant, ||9A' 11^00(^2-) :— ||9Ao|| ^00(^2) + ||9A'||ioo(^2') and ||(9fA'||^oo(^2-) :— 
II 5tA[, 11^00(^2) + liatA'lli, 00(^2). By similar arguments 

(4.9) ||J-*(A) - J*^*(A)|U2(i2) < Ch'/^ (||Ao|U.o(ia)||aA|U^(i2) + P|Uoo(i6)||Mo|Uoo(i2)) . 

Combining the above estimates together with the stability of J* we get 

(4.10) Kt < Ch'/^ (IIA'IUoo(^i) + ||atA'|Uo.(i2)) . 
Remark 6. If A is smooth, then 

(4.11) Kt < Ch (IIA'IUoc(Hi) + ||atA'|U,^(L2)) . 



From equations (2.17) and (2.4), we estimate the term Ks as 
(4.12) 



K.. 



f ({T'{A),dAA') ~ {rT'{A),rdAA'] 
Jm ^ 



< 2|| J-^(A)|U2(i2)||d^A' - J^d^A'|U2(i2) + 2\\T'iA) - J^J-'^(A)|U2(i2)|| J^d^A'|U2(i2). 
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Remark 7. By similar arguments as we used to bound K^ we get 

(4.13) K, < Ch'/^ (P'IUoo(^i) + WdtA'WL^^L-)) . 
If A is smooth, then 

(4.14) Ks < Ch (IIA'IUoo(^i) + \\dtA'U^^i^2^) . 

Summing up we get 
(4.15) 

implying consistency between ^[A] and 5"^ [A]. 

4.2. Consistency between 5^ [A] and 5'' [A]. The interpolated FEM action is given in equation 
(2.171, and the differential of it is 

DS''[A] ■ A' = DS-^[A]t ■ A' + DS-^[A]s ■ A', 

DS''[A]t • A' = 3fi V Mf^f>tr{j} (J-*(A) J*,(dAA')^ + h.c), 
(4.16) /../,' ' 

DS-^[A]s ■ A' = ^Y, Mf f, triJf {J'" (A)) Jf, {dAA' f +h.c.). 
fj' 

Let /( = {v, jViir+At, *r+At} be a temporal face, oriented from v — > jV — > jr+At -^ ir+At, 
and write ei = {(irjir}, 62 = {jrJr+At}; 63 = {jr+At, ir+At}, and 64 == {ir+At,«r}- For such a 
face, the constants Cgje = ~Ceet in the definition of J* take the value Ceie2 — C'ejes = Cesa — 
Ce4ei = I- Hence 

(4.17) 

4(dAA') - A',^ + A',,,^+K^ + A^,,^ + ^(K^,^o,eJ + [Ae,,yl^,,J) + 

^([^U'^ea] + [Ao,e,,A',3]) + -(K3, Ao,eJ + ^3,^0,64]) + 4([Ae4'^ei] + [Ao,e4,<])- 

In addition, we know that A^-^ + A^.^ = {dtA)f^ and v4o,e2 + ^o,e4 = {dAo)f^, implying that 
Aej = — ^ei + (c^t^)/t and Ao.ci = — Ao,e2 + (c'^o)/*, with similar relations for A'^^ and A'q ^^. 
This gives 

.^ ^g. 4,(dAA') = < + A;,^e2 + < + ^0,e4 + (K^^O.eJ + Ki,A[,_,J) + 

+ 0(AodtA' + ^dAo + A f^ A')/,, 

where e.g. O(v4odtv4')/t means 0((Ao)e,(dtv4')/J, e* S ft- 

Let / = {i, J, fc} be a spatial face, oriented from i ^ j ^ k, and write ei = {*, j}, 62 = {j, k}, 
63 = {^,i}- For such a face, the constants Cee' — —Ce'e in the definition of J^ take the value 

£162 ^^esei ^-^£263 g- -Li»:;ii»-»:^ 

J|(dA^') = A;^ + A',^ + A',^ + i(K,, Ae2] + [Ae4, A^J) + 

(4.19) 6 

-([<,Ae3] + [Ae4,<])+ -(K„^e3] + [^e2,<])- 

In addition, we know that A^.-^ +Aq.^ + A^^ = {dA)f, which means that A^^ = —A^-^ ^ A^^ + (dA)^, 
with a similar relation for A'^ , implying that 

(4.20) rjidAA!) = A,^ + < + < + ^(K,, Ae2] + [Ae4, A^J) + 0{AdA' + A'dA);. 
Again, ©(AdA')/ = 0{A^{dA')f), e e /. 
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Furthermore we need to calculate the differential of S^ [A] . The intermediate action is given in 
equations (2.281 and (2.291, and the differential of it is 



(4.21) 



DS^[A] ■ A' = DS'[A]t ■ A' + DS^[A]s ■ A', 
DS' [A]t • A' = 5R ^ Mf^fdi - ^* 



ftJl 



dr 



(F;,(A)-l)^+h.c. , 



/ dF^(A + tA') 
DS^[A]s ■ A' = SR^ Af^-/,tr f ^ _ (Ff (^) - 1)^ + h.c 



/,/' 



Thus we need to calculate 



(4.22) 



dF)^(A + rA') 



dr 



dF^JA + rA') 



r=0 



dr 



T=0 



We choose to calculate this by first writing the exponential functions in F as a single exponential 
using the BCH formula, and then using the formula for the differential, proposition [3] 

Let ft be as above. Then Fj (A) = UeJJe^UeJJei, where Ue^ = exp(^ei) and with similar 
expressions for Ue^, Ue^ and Ue^- By the BCH formula, proposition Bl Fl (A + rA') can be written 
as 



(4.23) 



F) (A + tA') ^ e'^'=i'^^^e^''''=2(^)e'^'=3^^)e'^''''=4(^) ^ (,^fMoi'Ao,^^,A,3,Ao..i)(r) ^ 



where A(t) = A + tA', and W/j(Aej,74o,e2 7^e3, ^o,e4)(''') is given from the recursion formula in 
proposition B] If we write Wf^{T) :— W^/j(ylei, ^o.ea: ^63,^0, 64 )(''') = X^^i^"' then the two first 
terms are 

di = Ae, (t) + Ao,e, (r) + Ae3 (r) + Ao^e4 {t) , 

(4.24) ^^2 = -Ki(T),Ao,e,(r)] + -[A,,{t),A,,{t)] + -[A,,{T),Aa^,,{T)] + 

We note from proposition Wl that dn oc C'(^A(r)"). If we again use A^^ — —A^,^ + {dtA)f^ and 
^0,64 = "^0,62 + (c^^o)/ti with similar relations for A' and A'q g^, then we get 



(4.25) 



Wf, (r = 0) = 4^ (-F*(A)) + ©(A^ + AdA)/, , 
dM//,(r) 



rfr 



T = 



J}^ (dAA') + ©(A^A' + AdA' + A'dA)^:^ , 



which again implies that 

(4.26) Fl (A) = e^f^ (^=°' = 1 + 4, [J"' (A)) + ©(A''' + AdA)/, . 

Using proposition |3] we get 

dFl{A + TA') 



(4.27) 

where 
(4.28) 



dr 



1-e 



= n(A) 



g-'^'i«',,(0) dW}^{T) 



ad 



Wf^iO) 



dr 



-adwjio) 



ad 



Wj{0) 



L ±0 

1 - -adwi(o) + ^arfvi/^.(o) 



+ 



and where we recall that adxY = [X, Y]. 
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Finally, combining the above estimates gives 
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(4.29) 

(F/,(A)-l)^i^/,(A) 



l-e""""'/*<"' dWf.ir) 



ad 



Wf, (0) 



dr 



4(^*(A))-'^ 



= Jj-,{T\A))" JJ^idAA') + 0{A'^A' + (A^ + AdA)dA' + (A^dA + (dA)2)A')/,/^'. 



T=0 

t t-rtfk\\H jt 1^ A'^ 1 /o/a4a/ 



Let / be as above. Then FUA) — UeiUe^Ue-i, where Uei ~ exp(ylej and with similar expres- 
sions for Ue2 and Ue^- By the BCH formula, proposition lil Ff{A) can be written as 



(4.30) 



F'JA + tA') = (.A,^(r) ^A,^(t) ^A^^(t) ^ gWj(Aei,A,,,A,3)(T)^ 



where A{t) — A+tA' , and Wf{Ae-^, A^^, Ae^){T) is given from the recursion formula in proposition 
H] If we write W^/(t) := Wf{Ae-^ , Ae2,Ae^){T) = X]^i dn, then the first two terms are 



(4.31) 



di ^A^^{t) + A^^ir) + A^^^{t), 



and we note that d„ oc 0{-A{t)'^). If we again use that A^^ = —A^^ — A^^ + {dA)f. with a similar 
relation for A^^, then we get 



(4.32) 



Ws(t = 0) = .//(^'^(A)) + 0(A3 + AdA)f, 
dWfit) 



dt 



t=o 



JfidAA') + 0{A^A' + AdA' + A'dA)f, 



which again implies that 



(4.33) FJ{A) = e^^(*=°) = 1 + J|(J"'(A)) + 0{A^ + AdA)f. 

Using proposition |3] we get 



(4.34) 



dFf{A + TA') 



dr 



T = 



FfiA) 



-ad 



g-aa„^(0) dWfir) 



ad 



Wf{0) 



dr 



T=0 



and by combining the above estimates gives 



{F^,{A)-irFf{A) 



-ad 



l_p-aa,vjm dWfir) 



ad 



Wf(0) 



(4.35) 



jpiF^iA)) 



HdWfir) 



dr 



dr 



T = 



T = 



0(AM')// 



:Jf(J-^(A))^J|(d^A') 

+ 0{A'^A' + [A^ + AdA)dA' + {A^dA + {dAf)A')ff, 
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Summing up, we get 
(4.36) 



hJ[ 



j_^-ad„,^^(„) ^yy^.^(^) 



ad 



Wf, (0) 



dr 



T=0 



h.c. 



= ^Y. ^'^Mi^' (4(J-*(A))^4(dAA') + h.c.) + 



ft,n 

+ ^Y1 ^^ftfi^^ (^ (^^^' + (^^ + A<dA)dA' + (A^dA + (dA)2)A')^^^, 
ftJi 
DS-^[A]fA'+ 

+ ^Y1 ^'^ftfi^^ (^ (^*^' + (^^ + AdA)dA' + (A\iA + (dA)2)A'^ 



/t./t' 



'ftfl 



and 



DS'[A], ■ A' ^ 5R^ A%,tr (Fj^iA) ifFjiA) 



fj' 



-ad 



l-e""""/'"' dWfir) 



ad 



Wf{0) 



dr 



h.c. 



= SR^ A%-tr (Jf (^^(A))^ j;(d^A') + h.c.) + 



(4.37) 



+ 5R 



^ ^f//'tr (o {A'^A' + {A^ + AdA)dA' + [A^dA + {dAf)A')^^\ 



fJ' 



DS-'[A]s-A'+ 

+ 3?^ Af//'tr (o {A'^A' + {A^ + AdA)dA' + {A^dA + {dAf)A'] 
fJ' 



ff 



If we assume that the mesh satisfies a CFL condition, i.e. there exists a constant C such that 

(4.38) 

then we can deduce the bounds 



< ^ < ^ < C, Vn, 



(4.39) 



\A,\ < Ch'/^AU^^L^-,, |Ao^,J < C/ii/2|lAo|U..(i.), 
\idA)f\ < Ch'/^dAh^^L.^, \{dtA)f,\ < C/.i/2||aA|Uo.(i.), 
\{dAo)f,\<Ch'^^dAoh^^L--), 



and we can conclude 
(4.40) 



\DS-'[A]-A'-DS'[A]-A'\<C{\\A'\\L^^H.) + \\dtA'\\L^^L-^)) 
implying consistency between 5'^ [A] and S^ [A] . 



/ii/2 Aei°°(iJi),atAei°°(L2 



A smooth, 



4.3. Consistency between S^[A] and ^'■'^[A]. The only difl^erence between S^ and S^ is the 
parallel transport operators in S^ introduced to make S^ gauge invariant. They are given as 



(4.41) 



U^^, = exp(Ay.y,) = 1 + Aj.^., + 0{A^), 

^0,/(r)/(r') = ^^P(\j{T)f{T')) = 1 + \f(T)f(T') + C(^o)> 



and the differential of these are proportional to A' and A'q respectively. Hence, by similar con- 
siderations as in the previous section we get exactly the same estimate as in equation (4.401, 
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i.e. 
(4.42) 

\DS^[A]-A'-DS^[A]-A'\ < C {\\A'\\,^^h^^ + HA'h^iL^)) If' ^ ^ ^°"^f'^'^'^ ^ ^°°^^'^ 

I ft A smooth. 

We summarize the results in a theorem: 



Theorem 8. Assume M is a bounded domain in M.^^^ . Then the SGT action (2.341, is consistent 



with the continuous Yang-Mills action (2.4|, with respect to the norm 

(4.43) ||A|| := ||A|Uoo(ffi) + \\dtA\\L^^L^), 

under the assumption that the above-mentioned CFL condition holds. 

As a consequence of the above estimates, we get the following estimate for the deviation of the 
SGT action S^ from the continuous action S, 



{AAA) \S{A)-S^{A)\< 




, A e L°°{H^),dtA e L°°(i2) 
, A smooth. 



5. Numerical convergence tests 

The preceeding sections have defined and proven consistency of the SGT action. However, due 
to the complexity involved in these quantities, we would like to include some numerical convergence 
tests as well. 

In our computer calculations, we focus on pure gauge theory with Q — SU{2), and used a 
four-dimensional cubic euclidean domain [0, 1]"* C K* with periodic boundary conditions. The 
lattice structure consisted of a three-dimensional simplicial mesh replicated at each discrete time 
value. The three-dimensional simplicial mesh consisted of a homogeneous arrangement of A'''^ 
identical cubic building blocks, each building block containing six tetrahedra as shown in figure 
|5.1[ Each such spatial mesh was replicated N times in the time direction to fill the four-dimensional 
domain, in accordance with the construction detailed in the previous sections. The lattice constant 
h is defined as the side length of each cubic building block, and also coincides with the time 
discretization interval. In the interest of simplicity we enforced temporal gauge, in which the 
temporal link matrices reduce to indentity matrices. 

The SGT action employs parallel transport matrices in order for gauge invariance to be re- 
spected. By defining the distinguished points of all spatial and temporal faces to coincide for as 
many pairs of faces as possible, we only need the parallel transport matrices for terms in the action 
involving pairs of temporal faces with no common nodes. 

In order to test convergence of the euclidean SGT, we compared the discrete and continuum 
action for several different choices of gauge fields for which the continuum action S |2.4| can be 
calculated exactly. We chose the following cases: 

(1) Gauge field oriented towards the x-direction in space and towards the generator t^ — icr^/2 
within su(2), with a sinusoidal time dependence, where a^ is a Pauli matrix. The only 
nonzero component of the gauge field A is 

Al{t,x,y,z):^—sm{2TTt), S = 1. 

(2) Gauge field oriented towards the y-direction in space and t^ within su(2), with a sinusoidal 
x-dependence. The nonzero component of the gauge field in this case was 

Al{t,x,y,z) := —sm{2nx), 5 = 1. 
" Ztt 

(3) A case with two nonzero components, 

Al := — sm{2TTy), Al := — sin(27ra;), S 



27r ' "" y ' 2tt ^ " 2 8(27r)^ 
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Figure 5.1. Elementary 3d mesh building block containing six tetrahedra, all of 
which share the single interior diagonal. 



1. 






1, 



5"= -. 



(4) A constant field that only contributes to the nonlinear term in the field strength, 

1 

2' 

The first case is insensitive to the spatial face mass matrix elements, while the second is insensitive 
to the spatial edge mass matrix elements which are used in the definition of the temporal mass 
matrix elements. In the first two cases, the nonlinear contribution to the continuum field strength 
F vanishes, and the action can be calculated analytically to be unity. In the third case, the 
nonlinear term survives, and the exact value of the action is 

1 1 

2 



S 



8(27r)4- 
In all cases, we measure the relative error of the discretized action, versus the lattice size, for 



lattices sizes N from iV = 4 to A^ = 32. The results are displayed graphically in figure 5.2 



We plotted the absolute value of the relative deviation of the numerical action from the con- 
tinuum action, as a function of the number of lattice sites per dimension. In all cases, the errors 
approach zero as the lattice resolution grows, so the action converges towards the continuum result. 
Using least squares second order polynomial fit, we determined that the relative error depends in 
the following way on the lattice constant h, 

Relative error w Ch?' ^ 

where C is some constant depending on the choice of gauge field. This is in accordance with the 
estimate 14.441 

Quantum SU{2) gauge theory Monte Carlo computer simulations using the SGT action will be 
published in a companion article |[8j| . 



6. Conclusions 

We have proposed a general formulation of lattice gauge theory on simplicial lattices. For any 
simplicial lattice of arbitrary shape, this action can by used for lattice gauge theory simulations or 
study of the classical equations of motion. Traditionally, lattice QCD simulations within physics 
have used a homogeneous mesh. Mesh refinement is a well established concept within the subject 
of FEM. We feel that it is well worth the effort to invenstigate the possibility of mesh refinement 
within classical and quantum gauge theory. Quantum lattice gauge theory simulations are very 
computer intensive. Therefore mesh refinement could be beneficial in cases where it makes sense 
to focus more computational effort on some subset of the simulation domain. 
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Figure 5.2. The relative error of the action versus the number of lattice sites 
per side N, for the actions 1, 2, 3, 4 described in section |5] The squares are the 
simulation data points and the solid lines are the second order polynomial fits. 
Errors are proportional to h^ in all cases. 

We have shown the consistency of this SGT numerical approximation to the continuous action, 
in the sense of approximation theory. The lattice gauge theory formalism is of such a complexity, 
that it makes sense to complement this theoretical proof with numerical "evidence". We have 
provided this for a few different cases of gauge fields, for which the action was shown to converge 
towards the continuum result as the grid fineness increased. 
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